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Is the Sullivan Process Compatible with QCD Chiral Dynamics? 
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We calculate the leading non-analytic quark-mass dependence in the moments of isovector quark 
distributions using heavy-baryon chiral perturbation theory. The results differ from what has been 
obtained from the Sullivan process in which hard scattering occurs through the virtual pion cloud 
of the nucleon. Our results provide useful guidance in formulating meson-cloud models consistent 
with chiral dynamics and can be used to constrain the extrapolations of the existing lattice QCD 
results to the physical quark masses 



The nucleon (proton and neutron) has a quark sea. 
For a long time, the leading picture for this sea comes 
from the meson-cloud model in which the nucleon can 
virtually dissociate into mesons such as pions, kaons, 
etc., plus a baryon core. A concrete realization of this 
model can be found in lepton-nucleon deep-inelastic scat- 
tering in which the (Bjorken-like) virtual photon scatters 
deep-inelastically off the meson cloud. The Feynman di- 
agram for this process is shown in Fig. 1 and is called 
the Sullivan process Many interesting results from 
the meson-cloud model have been obtained, which have 
shed important light on the dynamics of the quark sea 




FIG. 1. The Sullivan process in which the virtual photon 
scatters off the meson cloud in the nucleon. 

A crucial question about the meson-cloud model and 
the associated Sullivan process is to what extent the 
model can be justified from the underlying quantum 
chromodynamics (QCD) and its low-energy effective the- 
ories? Recently, A. Thomas et. al. pointed out the best 
place to check is the chiral properties of the model predic- 
tions They worked out, for the first time, the detailed 
leading non-analytic quark (pion) mass dependence of 
the parton distributions probed in the Sullivan process. 
This letter seeks to answer the above question by calcu- 
lating the same type of chiral logarithms using the low- 
energy effective theory of QCD: the heavy-baryon chiral 
perturbation theory [Q. As it turns out, there are subtle 
differences between the traditional approach in calculat- 



ing the Sullivan process and chiral perturbation theory. 
Our results can be used to improve upon the formula- 
tion of the Sullivan process and to study the quark mass 
dependence of the lattice QCD calculations. 

We are interested in the moments of the quark distri- 
butions in the proton, 



dxx n ~ l (q{x) + {-l) n q{x)) , 
dxx^ 1 (A.? (x) + (-If -1 Aq (xj) , 



dxx^'iSq (x) + (-l) n Sq(x)) 
(n = 1,2,3,...) 



(1) 



where q (q) is the quark (antiquark) spin-averaged dis- 
tribution, Aq (Aq) the hclicity distribution, and Sq (5q) 
the transversity distribution || . The variable (Feynman) 
x is the momentum fraction of the proton carried by a 
quark in the infinite momentum frame, and for simplic- 
ity we have suppressed the renormalization scale depen- 
dence. In QCD, the moments are related to the forward 
matrix elements of the twist-two operators 



Q 



qj^ 1 j 5 iD' 12 ■ --iD^q , 

through the relations, 
(PS \0^-^ \PS) = 2 (x n - 1 ) q P^ 1 ■ ■ ■ P^ ] , 
(PS O^-^ PSj =2{x n - 1 ) Aq MS (pi P^ 2 •••P A1 
PS O^ 1 -"" PS\ =2{x n - 1 ) Sq 



(2) 
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where (• • •) and [• • •] denote, respectively, the sym- 
metrization and antisymmetrization of the indices in be- 
tween. jPS 1 ) is the ground state of the nucleon with four- 
momentum and polarization vector S 7 ' (S 2 = — 1). 
M is the nucleon mass. All tensors are trace free. 
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To calculate the leading non-analytic quark-mass de- 
pendence of these moments in heavy-baryon chiral per- 
turbation theory 0, we need to construct a chiral ex- 
pansion of the quark operators in terms of the hadronic 
(pion and nucleon) operators with the identical symme- 
try properties. In this study, we focus on the isovec- 
tor combinations of the quark densities, up minus down 
quarks. To leading order in chiral power counting, we 
have, 
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where C n , C n , and C n are the corresponding moments 
in the chiral limit, gA is the isovector axial charge in the 
same limit, and the renormalization scale \i can be taken 
to be 471"/^ ~ 1 GeV. The relative coefficients between 
the results in the chiral limit and the leading chiral- 
logarithms are seen to be fixed. For (l^u-Ad' wn i cn 
is the physical isovector axial charge, our leading non- 
analytic correction agrees with the previous calculation 
M. {l) u _ d is the isospin charge which is protected from 
chiral corrections by the isospin symmetry. The next im- 
portant chiral correction is of type m^/(47r/ 7r ) 2 which is 
linear and analytical in quark masses. These corrections 
cannot be computed without additional parameters. 

The leading non-analytical quark mass dependence of 
(x n ) u _ d was calculated using the pion-cloud model (the 
Sullivan process) in ref. j|. Instead of the pre- factor 
(3g A + l) in front of the chiral logarithm, they found 
+7iV5 [a S' (Alll u Al2 • • • (u^t 3 vJ - ut 3 u) N Ag\. Although in the g A = 1 limit the two results agree, 

they are different in general. In particular, in the limit 
of a large number of colors (N c ), gA ~ N c [||. 

We find that the above discrepancy comes from the 
differences between the linear and non-linear formula- 
tions of chiral expansion. While the non-linear formula- 
tion used in our calculation has a simple and clear power 
counting scheme, the linear formulation does not. The 
Sullivan process is traditionally done in the linear sigma 
model with the following lagrangian 



N V (^ . . . y^n) ( UT; , u t _|_ u t T3U ) N 

+alVS" (A11 ^ 2 • • • v» n) (u f r 3 u - UT 3 ti t ) N , 
NS^v^ 2 ■ ■ ■ (ur 3 M t + v?t 3 u) N 
+0Nv (tl1 ■ ■ ■ (u^t 3 u - ut 3 u^) N , (4) 
NS [a v { >* l] ■ ■ ■ (ut r3il t + ut 3 u) N 



where N and N are the nucleon fields, w M is the nucleon 
four-velocity, and u — exp(j7r a r a /2/ w ) with pion fields 
7r a and decay constant f w = 93 MeV. a, (3 and 7 are 
unknown coefficients. For n = 1, there is also a lead- 
ing pion operator in 0^_ d , which has not been shown 
explicitly. 



*(« $ + gnNN((T + iT •7?7 5 ))* + 



(6) 



FIG. 2. Feynman diagrams contributing to the leading chi- 
ral logarithms in the nucleon matrix elements of twist-two 
operators. The dashed lines represent pions. 

The contributions to the leading chiral logarithms of 
the twist-two matrix elements come from three separate 
sources: first, the wave function renormalization of the 
nucleon states; second, the virtual-pion cloud as shown 
in the first diagram in Fig. 2; third, the pion tad-pole 
contribution as shown in the second digram in Fig. 2. A 
straightforward calculation yields, 
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n = l,2,--- 



(5) 



where represents the nucleon field in the linear repre- 
sentation and g-nNN is the nucleon-pion coupling. The 
pure meson sector of the lagrangian has been omitted. 
After spontaneous symmetry breaking, the above model 
in principle yields a Goldberg- Treimann relation g n NN — 
M/ fn with gA = l, In the literature however, g^NN 
in the linear-cr model result is usually replaced by the 
full Goldberg- Treimann relation g n NN = gA^/fw 

To make the linear er model result consistent with 
that of non-linear chiral expansion, one must add higher- 
dimensional operators to the lagrangian in Eq. @. As 
we have mentioned above, the interaction terms of dif- 
ferent mass dimensions contribute to the same order of 
chiral power is the main disadvantage of the linear for- 
mulation. In the present example, we need to add the 
following dimension-five term ]l0| , 

£' = g'ty ((t • 7?i $0 — <jt ■ i ■ 7r)75 

+f • (ff x i fin)) ^ . (7) 

To restore the full Goldberger-Treiman relation, g' ~ 
<7,4 — l. Setting gA = 1 in the result of Ref. (|] and 
adding the additional contributions generated from the 
above lagrangian, we find the chiral logarithms in the 
linear formulation coincide completely with the results 
above. 

Our results can be used to constrain the quark mass 
dependence of the matrix elements calculated in lattice 
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QCD Indeed the isovector combinations of the 

moments are easier to obtain since they do not require 
a computation of expansive "disconnected" diagrams. 
The results can also be used to improve the meson-cloud 
model to the point where it becomes consistent with the 
chiral physics of QCD. 

To summarize, we have calculated the leading non- 
analytic quark mass dependence of the twist-two, isovec- 
tor quark matrix elements in the nucleon. The results 
indicate that the traditional approach in calculating the 
Sullivan process needs to be modified to be fully consis- 
tent with the chiral dynamics of QCD. 

Note added in proof: After this work was completed, 
we learned a similar work by D. Arndt and M. J. Savage 
|^2| . For the same quantities computed, both papers 
agree. 
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